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ABSTRACT 



I use the fact that the radiation emitted by the accretion disk of supermas- 
sive black hole can heat up the surrounding gas in the protogalaxy to achieve 
hydrostatic equilibrium during the galaxy formation. The correlation between 
the black hole mass Mbh and velocity dispersion a thus naturally arises. The 
result generally agrees with empirical fittings from observational data, even with 
Mbh — 1O 6 M . This model provides a clear picture on how the properties of 
the galactic supermassive black holes are connected with the kinetic properties 
of the galactic bulges. 

Subject headings: Galaxies, galactic center, supermassive blackholes, velocity 
dispersion 



1. Introduction 



In the past decade, some observations have led to some tight relations between the 
central supermassive blackhole (SMBH) masses Mbh and velocity dispersions a in the bulges 
of galaxies. These relations can be summarized as log(Ms^/M Q ) = /31og(<r/200 km s _1 ) + a. 
For 1O 6 M < Mbh < 1O 9 M , the values of a and /3 ha ve been estimated seve ral times in the 



past 12 years: originally (a, ff) = (8.08±0.q 8, 3.75±0.3) flGebhardt et alJl2000h and (8.14±1.3 
4.80 ± 0.54) jFerrarese and MerrittJ lioooh . then (8.13 ± 06, 4 .02 ± 0.32) iTremaine et al. 
2002h. and more rece ntly (8.28 ± 0.05, 4.06 ± 0. 281 Hboosh. f8.12 ± 0.08, 4.24 ± 0.41) 
dGiiltekin etaDl2009h. f8.29 ± 0.06. 5.12 ± 0.36) jMcConnell et alJboilh and (8.13 ± 0.05, 
5.13 ± 0.34) (IGraham et al.l 120111 ) . Generally speaking, empirical fittings show that a ~ 8 
and (3 w 4 — 5. These relations correspond to all morphological type galaxies. One can 
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separate the fittings int o different groups such as the early- type and late-type or elliptical 
and spiral. For example, iMcConnell et all (boilh obtain (a, 0) = (8.38, 4.53) and (7.97, 4.58) 
for the early-type and late-type galaxies respectively if they are fitted separately. The slopes 
are shallower than the combined one (0 = 5.12). Moreover, the slope and the scatter of 
the Mbh — cr rela tion are still subject to debate, particularly at the low mass ends. Recently, 
Xiao et al.l (120111 ) obtain a new Mbh — & relation with low BH masses (below 2 x 1O 6 M ). 
They find a zero point a = 7.68 ± 0.0 8 and sl o pe = 3.32 ± 0.22, which indicate may 



be smaller for lower BH masses. Also, IWyithd (120061 ) obtained a better fit by using a log- 
quadratic form \og(M BH /M Q ) = a + /31og(cr/200 km s" 1 ) + 7'[log(a/200 km s" 1 )] 2 with 
a = 8.05 ± 0.06, = 4.2 ± 0.37 and 7' = 1.6 ± 1.3. Therefore, it is reasonable to doubt that 
the relation is not simply given by Mbh oc a" with just a constant for all galaxies. 



The Mbh — & relation has been derived by recent theoretical models ( 


Silk and Rees! 


1998|; 


Adams et al.12001 


MacMillan and Henriksen 


2002 


Robertson et al. 


2005 


Murrav et al. 


2005 


King 


2005 




McLaughlin et al. 


2006: 


Power et al. 


2011 




Navakshin et al. 


2012). However, 



these models contain various assumptions and fail to explain the relations in the small 
SMBH mass regime (0 ~ 3.3). In this article, I present a model to get an exact Mbh — & 
relation, which can explain the parameters a and in the empirical fitting in both small 
SMBH regime and apply to different types of galaxies. I use the fact that the strong radia- 
tion of the accretion disk of a SMBH can heat up the surrounding gas so that hydrostatic 
equilibrium of the latter is maintained. The cooling of the surrounding gas is mainly given by 
recombination, bremsstrahlung radiation and the adiabatic expansion of the gas. Without 
any other assumptions, the exact Mbh — & relation is naturally obtained. In the following, 
I will first present the details of the model. Then I will fit our model with the data of a and 
Mbh from 198 galaxies and show that it generally agrees with the empirical fitting. 



2. The Accretion model of supermassive black hole and the M B h — & relation 

It is commonly believed that all SMBHs accompany with accretion disks to emit high 
energy radiation during their formation. The luminosity of the disk is mainly come from the 
rest mass energy of the mass accretion. The accretion luminosity can be expressed as 

W = vMc 2 , (1) 

where 0.05 < r\ < 0.4 is the efficiency of the process. Let fsa — L disk /L Ed , where L Ed = 
1.5 x 10 38 (Mbh /M q ) erg s -1 is the Eddington limit of accretion, we have 

L disk = 1.5 x W™f Ed erg s" 1 . (2) 
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The accretion disk of SMBH provides a large number of photons to heat up the sur- 
rounding gas in the protogalaxy during the galaxy formation. Assume that the power is 
mainly transmitted to the protogalaxy within a scale radius R through radiation by comp- 
ton scattering and photoionization. The optical depth of the gas is r = na p hR <C 1, where n 
is the number density of the hot gas and a p h is the effective cross section of the interaction 
of photons and hot gas particles, which is closed to the Thomson cross section OThom for 
zero metallicity. Therefore, the total power that can be transmitted to the protogalaxy is 
just L disk T. In equilibrium, the heating rate is equal to the co oling rate by bre msstrahlung 
radiation Ab, recombination A#, and adiabatic expansion A a (Katz et al.lll996l ): 



L disk r = A B0 n 2 T l ^ 2 V + A R0 n 2 T - 3 



1 + 



T 



10 6 K 



-l 



dV 



(3) 



where A B0 = 1.4 x 10~ 27 erg cm 3 s -1 , A R0 = 3.5 x 10~ 26 erg cm 3 s -1 , T, p, V, are the 
temperature, pressure and volume of the gas within R, r espectively. The t e rm pdV/dt can 
be w ritten as pdV/dt ~ pV 2 ^{^kT '/ 'm g ) l l 2 and p = nkT ( Muno et al.ll2004t IChan and Chu 
20081 ). where 7 ~ 5/3 is the adiabatic index and m g is the mean mass of a gas parti cle. The 
Viria l relation of the effective total mass of hot gas M g and T within R is given by f lSarazin 
1988h 



kT = fi GM s m, 
J 3R 



(4) 



where f\ is the virial constant. After the galactic bulge is formed, assuming spherical sym- 
metry and by Virial theorem again, one can get 



a 



f GM g 



(5) 



where f% is another virial constant. From Eqs. (3), (4) and (5), and assuming nV = M g /m g , 
we get 

L disk = L lC r 3 + L 21 a 2S (l + L^o 2 )- 1 + L, 



where 



Ll 



AboA j 



1/2 



^21 



3^ 2 ml /2 fl /2 Gk^a ph ' 



(6) 
(7) 



-22 



fim g 



3 x 10 6 M 
2.1/ 3/2 

p5/2 



(9) 
(10) 
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Take /i ~ 1 (isothermal distribution) and /-j 
with Eq. (2), we have 



1/5 (ICappellari et al.l 120061 ) . and combine 



M, 



BH 



1O 8 M 



fm [6.2^00 + 6.3^(1 



10.4^oo) _1 + 0.37r"V 5 



200 J > 



'IV. 



where 0200 — cr/200 km s -1 . Therefore, assuming r ~ 0.005, the last term dominates for 
a > 60 km s -1 , which agrees with the observed range of /?. The fitting parameters fsd and r 
are mainly controlled by the empirical fitting parameters a and (3 of the ob servational data 
respec t ively. In Fig . 1, we get an empirical fi t by us ing the data obtained from I Greene and Ho 
( 120061 ); Kiao et al.l (1201 lh : iMcConnell et al.l (120111 ) . The effective cross section due to metal- 
licity may contribute to a factor of 2-3 in Eq. (11). By usi ng the cross sections of som e major 
metals (carbon, nitrogen, oxygen, silicon) calculated from iDaltabuit and Cox! (1 19721) and as - 
suming metallicity of a protogalaxy is about 10~ 3 solar metallicity (IJappsen et al.ll2009l ). 
the effective ccross section is 2 x 10~ 24 cm -2 , which is about 3<TThom- Also, the estimation 
of the Mbh is not too reliable for Mbh < 1O 6 M . Therefore, the fitting parameters are 
just an order of magnitude estimation. In Fig. 1, the functional form of Eq. (11) generally 
matches the observational data. The best fitted parameters are f^d = 50 and r = 0.005, 
with 4.9% rms error. By fitting with the form \og{MsH / M & ) = /3 log 0200 + a, we can get 
(a, (3) = (8.2,4.5) with 5.1% rms error. Therefore, two functional forms can fit the data 
equally well. However, if we neglect the first two terms in the right hand side of Eq. (11), 
the best fitted line is M##/1O 8 M = 4.5a 2 5 o with 12% rms error. As mentioned above, 
the first two terms are significant when Mbh or a is small. Therefore, the effective slope 
of the log Mbh — log a relation is shallower (f3 < 5). That means the effect of cooling by 
recombination and bremsstrahlung radiation should be considered especially in galaxies with 
low velocity dipersion. 

Our result is consistent with the recent observations which indicate that many supermas- 
sive black holes may in volve a long period of moderate super-Eddington accretion (/^ ~ 10) 
during their formation (IKawaguchi et al.ll2004l ; iBrian and Zhaoll2004l ; IWang et al.ll2008l ). On 
the other hand, t he central number density in the Milky Way is about 0.1 — 0.5 cm -1 
(IMuno et al.l 120041 ). which corresponds to r ~ 0.001 — 0.005 in the bulge. The best fitted r 
is also consistent with the observational data. 



3. Discussion 

In this article, I present a new model to explain the Mbh — & relation in galaxies. The 
Mbh — & relation is not simply a power-law form Mbh oc a 13 , but partially depends on 
cr 3 and a 5 . This exact form of the Mbh — & relation agrees with the recent observational 



-5 - 



data, especially in the small Mbh regime. In this model, we can obtain 3 < < 5 and 
a ~ 8, which is consistent with the observational data of the grouped galaxies (s mall SMBH: 



a ~ 7.7, & 3.3; early-typ e: a « 8.4, « 4.5; late-type: a « 8.0, /3 « 4.6) fjXiao et al. 



20111 : iMcConnell et al.ll201l[ ). In general, this model suggests that larger r and fEd result in 
smaller slope and normalization constant a in the relation. Thus, if similar galaxies have 
similar bulge structure and accretion disks, then the Mbh — & relation of this particular 
type may be tighter. It generally agree s with the observatio n that the Mbh ~ c relation in 



elliptical galaxies only is less scattered ( iGraham et al.ll201ll ). 



In this model, the evolution pattern of the supermassive black hole does not affect the 
function form of the relation. The only physics here is the energy balance of the gas between 
the heating by the radiation from accretion disk and the cooling by the free-free emission, 
recombination and the adiabatic expansion of the gas particles. If the black hole is still 
significantly accreting, the energy given out would be balanced by the cooling of gas, which 
gives the Eq. (11). When the black hole's activity is switched off, the relation between the 
kinematic properties of the bulge and the Mbh has already been established, which remains 
unchanged in Eq. (11). Therefore, the exact relation between Mbh and a can definitely 
apply in both active and non-active galaxies. 

All the parameters obtained (fEd ~ 10 and r ~ 0.001) are consistent with the theoretical 
estimation and observation. Generally, our result supports the moderate super-Eddington 
accretion during the SMBH formation. The variations of fEd and T within the groups of 
galaxies may result in the observed scatter in the Mbh~® fittings. All the above results arose 
from existing natural physical laws without any extra assumptions. This model provides a 
clear picture on how the properties of the galactic supermassive black holes are connected 
with the kinetic properties of the galactic bulges. 
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Fig. 1 . — Mbh,r versus a for 19 8 galaxies obtained from I Greene and Hoi (120061 ): Kiao et al. 
(boilh : iMcConnell et all (boilh . where M BH , 8 = M BH /10 8 M Q . The solid line is generated 
from Eq. (11) with f^a = 50 and r = 0.005. The dashed line is in the form log(Ms#/M ) = 
a + /31og(o"2oo) with a = 8.2 and (3 = 4.5. The dotted line is Mbh,s = 4.5cr| 00 . 



